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Abstract 



In this paper, we study the control theoretic properties of a couple of interacting spin 
I's driven by an electro-magnetic field. In particular, we assume that it is possible to 
! observe the expectation value of the total magnetization and we study controllability, 

observability and parameter identification of these systems. We give conditions for 
controllability and observability and characterize the classes of equivalent models which 



, have the same input-output behavior. The analysis is motivated by the recent interest 

I in three level systems in quantum information theory and quantum cryptography as 

■ well as by the problem of modeling molecular magnets as spin networks. 

^: 

1 Introduction 

^ ' In recent years, there have been several proposals to use three level systems, the so-called 

qutrits, in quantum information theory. The proposals concern the use of these systems as 
building blocks for protocols in quantum cryptography [12] and communication [7] as well 
^ , as for the encoding of two logic qubits [16]. They also have been used to study fundamental 

questions in quantum mechanics such as entanglement measures [5], [9], [14]. A study of 
control of three level systems was considered in [6]. From a quantum control perspective, a 
system of two coupled three level systems represents the next more difficult case after the 
well studied system of coupled spin ^'s [11], [19], [22]. Motivation to study these systems also 
comes from the problem of modeling molecular magnets. These novel materials [3], [4], [13], 
[20], [21] are of interest in many applications as nanosize magnets as well as for fundamental 
studies in quantum mechanics and biology. They are modeled as networks of interacting 
spins. Spin I's are a very common example of three level systems. Examples are the nuclear 
spins of the naturally occurring isotopes ^Li, "^H, ^^iV. 

We shall study the control-related properties, namely controllability, observability and 
parameter identifiability, for a pair of interacting spin I's particles. To be more specific, we 
shall consider an Heisenberg spin model with Hamiltonian given by 

H{t):=tiA + B.,u.,{t) + ByUy{t) + B,u,{t)), (1) 



1 



with 



-iJi2{T, 



j=x,y,z '^j 



(2) 

By := (71(7,, (g) 1 + 72I (S) (T^), for v = x, y, or z. 

Here J12 is the exchange constant, 71 and 72 are the gyromagnetic ratios of particle 1 and 2, 
respectively; Ux,y,z are the x, y, z time- varying components of the input electro- magnetic field; 
1 is the 3x3 identity matrix. ax,y,z are the spin matrices spanning the three dimensional 
representation of su{2) [23] 
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The total magnetization for the state p in the direction v — x,y, z is given by 



(3) 



(4) 



(5) 



(6) 



where Sj'~''^ = (Ty C?) 1 + 1 ® (Ti,, V = x,y,z. Recall that the density matrix (state) of the 
system p satisfies the Liouville's equation [23] 



p^-i[H{t),p]. 



(7) 



We are interested in the Heiscnberg Hamiltonian (1) because we have in mind applications 
to spin Hamiltonians modehng the dynamics of molecular magnets [3], [4], [13], [20], [21]. 
However, the methods presented in this paper can be generalized to different types of coupled 
three level systems as for example two spins I's with interaction different from the one 
modeled in (2) or cases where one component of the magnetic field is held constant. The 
main tool is a Cartan decomposition of the Lie algebra su{3), described in Section 2, which 
gives a decomposition of higher dimensional Lie algebras constructed with tensor products of 
matrices in su{3). We begin by stating the definitions concerning controllability, observability 
and parameter identification with reference to the system we want to study. 



Definition 1.1 An n— level quantum system is controllable if it is possible to drive the 
evolution operator to any value in the special unitary group SU (n) . 

ControUabihty can be checked [18] by verifying that the Lie Algebra generated by the ma- 
trices defining the dynamics (^4, B^, By, B^ in (1), (2)) contains su{n) (in this case su{9)). 
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Definition 1.2 Denote by p{t, u, po) the trajectory corresponding to an initial state po and 
control(s) u. Let S be the matrix corresponding to the output of the system Tr{Sp{t, u, po)) 
m our case V = X, y, z in (6)). Then the system is observable if Tr{Sp{t, u, po)) = 

Tr{Sp{t,u, p'q)), for every t and control u, implies po = p'q. 

This definition of observability refers to identification of the initial state by a measure of the 
expectation value of an observable. Observability for quantum mechanical systems in these 
terms was studied in [10]. If C is the dynamical Lie Algebra (generated by A, B^^y^z above 
for our system), an n— level system with output S (assumed w.l.g. traceless) {S^^^ in our 
case) is observable if and only if 

V :^®f^QadciS ^ su{n). (8) 

Here and in the following ® denotes the sum of vector spaces (not necessarily direct sum). 
In the controllable case su{n) C C and (8) is verified. 

Now consider two models E and E' of Heisenberg spin I's. These models may differ by the 
parameters J12 and 71^2- They may also have different initial states say po and pg. Therefore 
we often consider the pair (S, po) and the pair (S', pp). We investigate whether it is possible 
to distinguish state and parameters by an experiment involving control with an input field 
and a measurement of the output 5"™^. This problem is motivated by recent results on 
the isospectrality of Heisenberg Hamiltonians which showed the impossibility to distinguish 
the parameters in the Hamiltonian by a measure of thermodynamic properties [24] . We call 
p{t) := p{t,Ux,Uy,Uz, Po) a general trajectory for (S,po) and p'{t) := p{t,Ux,Uy,Uz, p'q) the 
corresponding trajectory (with the same control) for (E', p'q). We give the following definition 

[!]■ 

Definition 1.3 Two pairs (E,po) and (E',Pq) are equivalent if 

Tr(5™^p(t))=Tr(5™VW), (9) 

for every trajectory p and corresponding (with the same control) trajectory p'. 

The question of whether or not it is possible to distinguish two models using a reading of the 
total magnetization will be posed by describing the classes of equivalent pairs model-initial 
state. If Po and p'q are scalar matrices so are p{t) and p'{t) for every t so the outputs (9) are 
identically zero independently of the model. We shall exclude this degenerate case in the 
following treatment. 

The rest of the paper is organized as follows. In the next section, we describe some 
properties and a decomposition of the Lie algebra su{3) that will be used in the following. 
The question of controllability and observability is tackled in Section 3, where we prove that 
the system is controllable and observable if and only if 71 7^ 72 and J12 7^ 0. In Section 4 
we give a description of the classes of equivalent pairs which, we prove, consist of only two 
elements. Some concluding remarks are presented in Section 5. 
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2 Properties of the Lie algebra su{3) 



The Lie algebra su{3) appears in several areas of quantum physics. As a result, it has 
been extensively studied in the mathematical physics literature (see e.g. [8], [15]). We 
describe here some properties that are important for our treatment. We consider a canonical, 
orthogonal basis of su{3) given by the three matrices 
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(compare with (3)- (5)), and the matrices 
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(11) 



(12) 



The matrices (Tx,y,z span a Lie algebra isomorphic to su{2) which we denote by S. In 
particular we have 



2a. 



[(7y, a. 



ar. 



a-,, 



(13) 



The matrices R,Q,T,V,U along with multiples of the 3x3 identity 1 span the orthogonal 
complement of <S in u(3), which we denote by S-^. The following tables summarize the 
remaining commutation relation for su{3) in terms of the basis elements we have defined 
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(14) 



Prom these tables, it follows 

[S, S] = S, [S^, S] = S^/span{il} [S^, S^] = S. 
We also have that for any matrix in the set {R,Q,T,V,U}, say L 

oo 

^adsL — S'^/span{il}. 

k=0 



(15) 



(16) 
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The anticommutation relations are summarized in the following tables 
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We have 



t{S,S} 



(18) 



In the following, we denote by o" a generic element of iS and by S a generic element of 
iS-^. Therefore, a and S are Hermitian matrices. The decomposition of ^(3) which we have 
introduced in this section has consequences for decompositions of higher dimensional spaces. 
We shall use this in the following sections, in particular in Section 4. 



3 Controllability and Observability 

The system of two interacting spin I's, if the gyromagnetic ratios are equal, has dynamical 
Lie algebra isomorphic to su{2) or -^(2) according to whether or not J12 is equal to zero. In 
the case J12 = we also have a Lie algebra isomorphic to su{2) even in the case of different 
7's. The only nontrivial case is when 71 7^ 72 and J12 7^ 0. In this case, we have the following 
Theorem. 

Theorem 1 // 71 7^ 72 and J12 7^ 0, the system is controllable namely the dynamical Lie 
algebra is equal to su(9). 

Proof. We have to prove that, by calculating (repeated) Lie brackets of the matrices A, B^^y^z, 
we can obtain all the matrices of the form iC ® D where C and D vary in the orthogonal 
basis of m(3) described in the previous section, except the 9x9 identity. By repeated Lie 
brackets of the B^^y^z and using a determinant of Vandermonde type of argument similar to 
the one in Lemma 4.1 of [2], we obtain all the matrices of the form ia®l and il ® a. Then, 
using the Lie bracket of these matrices with A several times, we obtain also all the matrices 
of the form iai ® a2- To obtain the other elements we proceed as follows: We calculate 
[ia, ® CTxjiaz (8) ay] (see Table (17)). This gives a multiple of il ® cr^ (which is already in 
the dynamical Lie algebra) plus a multiple of T (g) cr^, with T defined in (11). Prom this, 
taking the Lie brackets with elements of the type il<Sia and icr (g) 1, using (16) we obtain all 
the matrices of the form iS <^ a and analogously we can obtain all the matrices of the type 
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ia (g) 5". To conclude the proof of controllability we only have to prove that we can obtain all 
the matrices iS S except the 9x9 identity. Notice that since su{3) is a simple Lie algebra 
[sn(3), su{3)] = su(3). Therefore given C in su{3) we can choose two matrices M and N such 
that [M.N] = C. Using the well known fact (see e.g. [23]) that Y.j=x,y,z { ~ "i-^jY = 2x1, 
we calculate 

[M®-iaj,N ®-iaj\^ ^ C ® {-iajf ^2C ®1. (19) 

j=x,y,z j=x,y,z 

Analogously we can see that we can generate all the matrices 1 C with C e su{3). Now, 
since C is a general matrix in su{3) we can obtain all the elements of the type K ®Y with 
Y e su{3) (or K e su{3)) and K (or Y) in the orbit ®feLo^^L(3)^'^- However this orbit is 
equal to su{3) (it is a nonzero ideal in su{3) and therefore it must be su{3) itself since su{3) 
is simple). This concludes the proof. □ 

In the case 71 7^ 72 J12 7^ the system being controllable is also observable. In all 
the other cases, the space V defined in (8) is different from su{9). In these cases, initial 
density matrices which differ by a matrix in V"*" cannot be distinguished and the system is 
not observable. 



4 Parameter Identification 

We now characterize the classes of equivalent pairs model-initial state. In other terms, we in- 
vestigate what can be said concerning the parameters of the system by experiments involving 
control with an external electro-magnetic field and measurement of the total magnetization. 
We shall assume that we are in the controllable (and therefore observable) case, namely 
we know that 71 7^ 72 and J12 7^ 0. We state and prove the main result of this section 
in the following Theorem 2, where we characterize the classes of equivalent models. In the 
following we mark with a prime ' every symbol concerning system S'. We first give three 
prehminary results that can be proved as in the case of networks of spin |'s treated in [1]. 
For completeness we give self contained proofs and some additional considerations in the 
Appendix. 

Lemma 4.1 If, for every trajectory of E, p, and corresponding trajectory of E', p', we have 

Tr{Sp)^Tr{S'p'), v^x,y,z, (20) 

for some pair of matrices S and S', then for every F, F :— ads^^adBj^ ■ ■ • adsj^S, and 
corresponding F', F' :— ads', ads'. ■ ■ ■ adB'. S', (ji, jV £ {x, U, z} or Bj — A), we have 

Tr{Fp) = Tr{F'p'), (21) 
for every pair of trajectories p and p'. 
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Lemma 4.2 Let (S,po) and (S',Pq) be two equivalent models. Then up to a permutation 
of the indices 

71,2 = 7i,2, (22) 

and for every a & iS 

Tr{(T®lp{t))^Tr{(7®lp'{t)), Tr{l ^ ap{t)) ^ Tr{l ^ ap'{t)). (23) 

Lemma 4.3 Assume two models (S,po) and (S',Pq) are equivalent. For every pair of ma- 
trices S and S' such that 

Tr{Sp)^Tr{S'p'), (24) 

we also have 

Tr{[S,a^l]p)^Tr{[S',a^l]p') Tr{[S,l ^ a]p) ^ Tr{[S',l ^ a]p'). (25) 

We define now two orthogonal subspaces of isu{9): T which is spanned by elements of 
the type cxi ® (J2 and 5*1 ® S2 (namely the factors of the tensor product are both in iS or both 
in iS^), except the identity, and T-*- which is spanned by mixed type of elements namely 
elements of the type a <Si S and S <Si cr. We shall use this decomposition of isu{9) (which 
induces a decomposition of su{9)) in the following treatment. The induced decomposition 
of 5^(9) is a Cartan type (see e.g. [17]) of decomposition as stated in the following Lemma. 

Lemma 4.4 

su{9) = il®il^, (26) 

with 

\iI^,iI^]CiI^, (27) 

[iJ^,iJ]CiJ, (28) 

[iX,iJ]CiJ^. (29) 

Proof. To verify (27), we consider a Lie bracket ® Si,(T2 ® 5*2] and prove that it is 
orthogonal to elements of the form a-^ ® as well as to elements of the form S3 S4. To do 
this, we rewrite [cxi (S> 5*1 , (72 (S> 5*2] as 

[CTI ® 5*1, (T2 ® 'S'2] = 0"l(T2 (8) S1S2 — (72(71 (g) 5'2S'i. (30) 

We can decompose (7i(72 <S> S1S2 as 

(71(72 (8) ,5i,52 = ^([c7i,(72] + {a,, a 2}) ® ([^1,^2] + {^1,^2}). (31) 

From this expression, using (15) and (18), the only term that is not perpendicular to (73(8x74 
is [di, (72] ® [5*1, 5*2]. Doing the same thing for the second term on the right hand side of (30), 
one obtains that the only term which is not perpendicular to (73 <S> (Ta is [(72, cti] ® [^'2, 5*1]. 
But these two terms cancel. Analogously one proves orthogonality to matrices of the type 
<S'3 ® S4. To conclude the proof of (27) one has to prove orthogonality of terms of the 
form [(7i (S> 5*1, S'2 ® 0-2]- This is obtained using similar arguments. Also similar arguments, 
considering all the sub-cases, prove (28) and (29). □ 
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Remark 4.5 The argument in the above Lemma can be generahzed to deal with decompo- 
sitions of sm(3"), for every n > 1. One can define a subspace of isu{3"') of tensor products 
of matrices of the form c ® 5 ® ■ ■ ■ ® o" with an odd number of ci's and a complementary 
space with an even number of cr's. Call these subspaces Tg and Xg respectively. Then one 
can show, by induction on n that 

[iIo,iIo\ Qi^o, [iIo,i^e\ Qi^e, [ilci^e] Qi^o, (32) 

and 

C Xe, {Xo,Tg} C Xg, {Xe,Xg} C Xg. (33) 

In fact for n = 1 (32) and (33) follow immediately from (15) and (18). For n > 1, (32) 
follows by writing 

[A C L>] = ^{{A, C} (g) [B, D] + [A, C] {B, D}), (34) 

and applying the inductive assumption to all the factors in this expression and considering 
all the sub-cases. Analogously one can prove (33). 

The following is the main theorem of this section. 

Theorem 2 Two controllable pairs Model-Initial State E(n, J12, 71, 72, Po), S'(n', 7(2, 7i, 72, p'q) 
are equivalent if and only if (up to a permutation of the indices) 



1. 


71 = 7i 




2. 


72 = 72 




3. 


\Ju\ = 1 


V 1 
•^121 


4- 


IfJl2 = 


7(2 then Po = Po- If J12 = ^Ji2> denote by pi (p[) the component of p (p') in 



and p2 (p'2) the component of p (p' ) inX then pi(0) = p'i(O) and ^2(0) = —^2(0). 



Proof. Assume first that the two pairs are equivalent. We have, from Lemma 4.2, that, up 
to a permutation of the indices, 71,2 = 7i 2- 

Consider now the following procedure to generate a basis of su{9). Start with A, il® c 
and icT ® 1 at Step 0. At step n take the Lie brackets of the matrices obtained at step 
n — 1 with A, il ® a and ia ® 1. By controllability, the procedure generates a basis of 
sw(9). Moreover every element we calculate belongs to either iX or iX-^ and there are no 
combinations. This follows by induction on the step and applying Lemma 4.4. We can 
repeat the same procedure starting with A\ il ® a and ia ® 1. Let F and F' be two 
corresponding matrices obtained at a step d>l. We have F = J12F and F' = J12F for the 
same F and with k odd for F{F') e iX and even (not zero) for F{F') e iX-^. This is true 
at Step 1 and follows by induction for elements obtained at the following steps by applying 
Lemma 4.4. Now notice that elements obtained from Step 1 on also span all of su{^) as 
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well as the elements obtained including Step 0. This is because, if we call Ri the subspace 
spanned by elements obtained at Step 1, the elements from Step 1 on span the vector space 

(B^^k^^ ^_i^^^Qad'^^adij?^ ■ ■ ■ ad'^^Ri, where Ti,T2, ...,Tr are in the set {A,il ® a^ia ® 1}. It 

follows from an application of the Jacobi identity that this is equal to (BkLo^d^Ri, where C 
is the Lie algebra generated by {A, il ® a,ia ® 1} which by controllability is sm(9). So this 
is equal to ©fcLo'^'^L(9)-^i which is a nonzero ideal in su{^) and therefore sm(9) itself since 
sw(9) is a simple Lie algebra. The same argument holds with A' replacing A. In conclusion 
we have, by applying Lemmas 4.1, 4.2 and 4.3, for any F El 

J^^TriFp) = J'^Tr{Fp'), (35) 

with k odd and, for any F E 1^, 

J^,Tr{Fp) = J[^Tr{Fp'), (36) 

with k even. In particular, by applying (36) for F = 1 (g) cr and comparing with (23) of 
Lemma 4.2 we obtain | J12I = |</i2l- 

The proof goes now in an analogous way to the case of spin | treated in [1]. We have 

two cases: If J12 = J(2' we have A = A', B^^y^z = B'^yz- ^^^^ ^^^^ since the systems 
are observable and we have the same input-output behavior then we must have po = Po- H 
J12 = —Ji2i then from the above discussion we have 

Tr{Gp^)^Tr{Gp',), VG e X^, (37) 

and 

Tr(Gpo) = -Tr{Gp'^), e I, (38) 

so the components of po and p'q in are equal while the components in I are opposite. 

To prove the converse of the theorem, the only nontrivial case is when J12 = —J'u- 
this case let us write the equation for p as 

p=[A + B{t),p], (39) 

and the equation for p' as 

p'=[-A + B{t),p']. (40) 

We can write p (p') as p := pi + p2, p' := p'l + p'2 with pi(') G and p2(') G X. Using 
relations (27), (28), (29) of Lemma 4.4 and noticing that A e il while B{t) e iX-*-, for every 
t, we can write the differential equations for pi and p2 as 

p, = [B{t),p,] + [A,p2] (41) 

P2 = [A,pi] + [B{t),p2], 

and the differential equation for p[ and P2 ^.s 

p[^[B{t),p[] + [A,p',] (42) 
p'2^[A,p[] + [B(t),p',]. 
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Combining these equations we obtain a differential equation for pi — p'^ and P2 + P2- 
particular, we have 

pi - p\ = [B{t), pi - p[] + 2[A, p2 + p'2] (43) 
P2 + P2 = 2[A, pi - p'J + [B{t),p2 + p^]. 

From equations (43), it follows that if pi(0) = p'i(O) and P2(0) = — P2(0), then pi(t) = p'i(t) 
and P2(i) = — P2(^)> f*^^ every t, and for every control B(t). In particular, since Tr{Sj'^'^ p) = 
Tr[S^'~'^ Pi) and pi = p'^, the two models are equivalent. □ 

5 Conclusions 

We have presented a control theoretic analysis of a system of two coupled spin I's. In 
particular, this concerns the controllability, observability and identifiability properties of 
this model. A Cartan decomposition of the Lie algebra sit(3) induces a decomposition of the 
Lie algebra sm(9) which plays a fundamental role in the control theoretic properties of this 
system. A similar situation was found in [1] for general networks of spin | and it is likely to 
appear for other type of networks of spins not necessarily equal to | or 1. We believe that 
the methods of analysis developed in this paper can be generalized to include for example 
different forms of the interaction, models where one or more components of the input field are 
held constant, networks with more than two spins (cfr. Remark 4.5). We have proved that if 
(and only if) the gyromagnetic ratios are different and the coupling constant is not zero the 
system is controllable and observable. In this case, we have characterized the set of equivalent 
models that give the same input output behavior. Our results are motivated by the problem 
of identifying the unknown parameters in molecular magnets through experiments involving 
driving the system with an input field and measuring the total magnetization. The analysis 
is also instrumental to the design of controls which will be considered in further research. 
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Appendix: Proofs of Lemmas 4.1, 4.2, 4.3 

Proof of Lemma 4-1 

Proof The proof can be obtained by induction on the depth of F and F' defined as the 
number of the operations ad in its calculation. For depth (21) is the same as (20). Now 
assuming that F has depth — 1, we can write for every r and t (considering a trajectory 
corresponding to controls identically zero from a certain instant on) 

Tr(Fe^V(T)e-^*) =Tr(F'e^V(T)e-^'*), (44) 

which taking the derivative with respect to t at zero gives 

Tri[A,F]pir))^Tri[A',F']p'ir)). (45) 

Analogously one can obtain (with appropriate constant control) 

Tr([A + F]p(r)) = Tr([A' + F']p'(r)), (46) 

which combined with (45) gives 

Tr([B,,y,,, F]p(r)) = rr([S;^,„ F']p'(t)). (47) 

□ 

Lemma (4.1) can be generalized as follows. For E and S' we can construct a basis for 
the dynamical Lie algebra starting from A, B^^y^z or A', B'^ y^^ and at each step calculating 
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the Lie brackets of the elements obtained at the previous step by A, B^^y^z or A', B'^ y^^. 
Consider the depth of the element of the basis as the number of Lie brackets calculated. The 
generalization consists of noticing that if (20) holds for some S and S' it also holds for [L, S] 
and [L', S'] where L and L' are elements of the basis of the dynamical Lie Algebra obtained 
the same way just replacing the A, Bx,y,z with A', B'^ y^^. This is true for every element of 
depth from Lemma 4.1. Now assume it is true for elements L and L' of depth d—1. Prom 
the Jacobi identity we have 

Tr{[[B, L], S]p) + Tr([[L, S], B]p) + Tr([[5, B], L]p) = (48) 

Tr([[B', L'lSy) + Tr([[L', SlB']p') + Tr([[S', B'], L']p') = 0, 

for some B in the set A, B^^y^z and corresponding B' . Now the second terms of the two sides 
are equal by applying the inductive assumption and Lemma 4.1. The same thing is true for 
the third term where we apply first Lemma 4.1 to obtain Tr{[S^ B]p) = Tr{[S' , B']p') and 
then the inductive assumption on L (with S replaced by [S, B] ) . Therefore the first terms 
are also equal. This facilitates the proofs of Lemmas 4.2 and 4.3. 

Proof of Lemma 4 ■ 2 

Proof. By performing (repeated) Lie brackets of -B^, By and Bz, it is possible to obtain all 
the matrices of the form 'f'^ia ® 1 + 72^1 ® cr, with k = 1, 2, ... (cfr. Lemma 4.1. in [2]). The 
corresponding matrices for E' are Yi'^ia <8) 1 + 72 '^il (8) cr. Now, starting from 

rr(5™^p)=rr(5™V), (49) 

and taking the Lie bracket with the matrices above obtained, we have 

7i^Tr((T ® Ip) + 72rr(l ® ap) = i^Tr{a ® Ip') + 72''Tr(l ® ap'), A; = 0, 1, 2, ... (50) 

Since Tr{a ® Ip) is not zero for every trajectory p (unless p is a scalar matrix which is a 
case we exclude), the only possibility for (50) to be verified is that the determinant 



D 



7i 72 7i 72 

ii ii 1? a 

il il ii ii 



(51) 



is equal to zero. But this is a Vandermonde determinant, therefore we need two of the 7's 
and 7"s to be equal. Up to a permutation we can choose 71 = 7i. We can now use the same 
Vandermone determinant type of argument starting from 

-l\{Tr{a ® Ip) - Tr{a ® Ip')) + 72^^(1 ® (^p) - l2Tr{l ® ap') = 0, (52) 

to conclude that 72 = 72 and that (23) holds. □ 

Proof of Lemma 4-3 
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Proof. As in the previous Lemma we obtain 

^'[{Tr{p[S, a 1]) - Tr{p'[S', a^l])+ j^{Tr{p[S, 1 a]) - Tr{p'[S', 1 ® 1]) = 0, A; = 1, 2, ... 

(53) 

This since 71 7^ 72, we obtain (25). □ 
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